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THE SOLUTION OF THE LAMINAR -BOUNDARY-LAYER EQUATION FOR 
THE FLAT ELATE FOR VELOCITY AND TEMPERATURE FIELDS 
FOR VARIABLE PHYSICAL PROPERTIES AND FOR THE 
DIFFUSION FIELD AT HIGH CONCENTRATION* 

By H. Schuh, 


SUMMARY 


In connection with. Pohlhausen’s solution for the temperature field 
on the flat plate, a series of formulas were indicated "by means of 
which the velocity and temperature field for variable physical charac- 
teristics can he computed hy an integral equation and an iteration 
method based on it. With it, the following cases were solved: On the 

assumption that the viscosity simply varies with the temperature while 
the other fluid properties remain constant, the velocity and tempera- 
ture field on the heated and cooled plate, respectively, was computed 
at the Prandtl numbers 12-5 and 100 (viBcous fluids). A closer Btudy 
of these two cases resulted in general relations: The calculations 

for a gas of Pr number 0 .7 (air) were conducted on the assumption that 
all fluid properties vary w^th the temperature, and the velocities are 
low enough for the heat of friction to.be discounted. The result was 
a thickening of the boundary layers, but no appreciable modification 
in shearing stress or heat-transfer coefficient. The effects of 
density and viscosity or density and heat conductivity have opposite 
effect for velocity and temperature field and almost cancel one 
another. Formulas allowing for the heat produced by the friction were 
indicated, but no calculations were carried through in view of the 
already existing report by Crocco . The methods of solution developed 
here were finally applied also to the case of diffusion of admixtures, 
where at higher concentration finite transverse velocities occur at the 
wall. 


*"Uber die Losung der laminaren Grenzschichtgleichung an der 
ebenen Platte fur Geschwindigkeits- und Tempera turf eld bei 
veranderlichen Stoffwerten und fur das Diffusionsfeld bei hoheren 
Konzentrationen. Zentrale fur wissenschaftlich.es Berichtswesen der 
Luftfahrtforschung des Generalluftzeugmeisters (ZWB) Berlin-Adlershof, 
Forschungsbericht Nr. 1980, August 18, 1944. 



2 


NACA TM 1275 


I. INTRODUCTION 


The laminar-boundary-layer equation for the flat plate lined up 
with the flow and with constant fluid properties was solved by Prandtl 
(reference 1) and Blasius (reference 2) for the flow field and by 
Fohlhausen (reference 3) for the temperature field. Since temperature 
and velocity field coincide when kinematic viscosity (V) and tempera- 
ture conductivity (a) are identically equal ^Pr = ^ = 1 ^ , Pohlhausen 7 s 

formula for the temperature field contains a solution for the velocity 
field also in the form of an integral equation. Piercy and Preston 
(reference 4), proceeding from a rough approximation, indicated that, 
with the aid of this integral equation and an iteration method, the 
well-known BlaBius solution can he obtained in a few steps. This 
method of solution has the advantage of being simple and. ^requiring 
relatively little time. It is - as is shown in the following - 
particularly suitable for boundary-layer calculations involving 
variable fluid properties, because a first, and usually fairly close, 
approximation, is already available in the solution for constant fluid 
properties . * 

Crocco (reference 5) and von Kdrman and Tsien (reference 6) (see 
also reference 9, 10) computed velocity and temperature field for 
variable fluid properties. In both reports, the differential equations 
are put in a different form from the elsewhere conventional boundary- 
layer calculation by changing to new variables. Crocco obtains two 
simultaneous differential equations of the second order which he solves 
for a gas with the Prandtl number Pr = 0.725 (air) . Ton hinndn and 
Tsien treat the case of Pr = 1 and have to solve only one differential 
equation, since then the temperature is related in a simple manner to 
the velocity. 

In the following, it is shown that a number of boundary-layer 
problems for the flat plate pan be solved in a comparatively simple 
manner, involving merely quadrature, by means of the cited integral 
equation and an iteration method. 


Ill SOLUTION OF BOUNDARY-LAYER EQUATION FOR VARIABLE PEYSICAL PROPERTIES 


the 


The boundary-layer equations for velocity and temperature field at 
flat plate at variable density read (reference 3) 


chi 

pU ^ + 




( 1 ) 
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+ Sigh . 0 ( la ) 

01 dy 

v“i + v T | - 4( x | ) (2) 

a 

with, u and v the speed in flow direction and at right angle to It, 
T, the temperature, 1 the distance from the plate leading edge, 
y the distance from the wall, p the density, u the viscosity, Cp the 

specific heat, and \ the heat conductivity. In the equation for the 
temperature field, the heat produced "by friction is, at first, not 
taken into account; as long as the speeds are not excessive and the 
temperature differences not too’ small, this Is justified. 


For constant density, equations (1) and (2) can he reduced to an 
ordinary differential equation (reference 3) on the assumption that u 
and T are a function merely of the one (dimensionless) 


Since the density depends only on the tempera- 


coordinate | = 

ture, the Idea suggests itself that the same simplification Is possible 
also for variable density. We put 


m(|) 


T - T 0 
Tq - T 0 


6(1) 


1 ■"H'S 


o) 


where TJ is the velocity at the edge of the boundary layer, 

T 0 and T-|_ the wall temperature and the temperature at the edge of the 

boundary layer, respectively. The quantity v jj. in the dimensionless | 
denotes the kinematic viscosity for the fixed temperature T^., for 

which in suitable manner the wall temperature (k = 0), or the tempera- 
ture at the edge of the boundary layer (k = 1), Is chosen. The 
boundary conditions for flow and temperature field read 


II 
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where the subscript k denotes the density at temperature T^., gives 
by (la) 

M ( puS 'll ou iS ) (S) 

hence, by (l), after introduction of (5) and (6), the more suitable 
form 



From (7), regarded as differential equation for the quantity cp^- and 

f temporarily as a known function of g, the following expression for 
oi is derived 



This disposes of the integration constant from consideration of the 
boundary condition (4) . Likewise, there is afforded for the dimen- 
sionless temperature 0 the expression 4 — - 



where Pr^ is the Prandtl number with the density at temperature T^.. 

For constant density (tp ■ y = X = 1), velocity and temperature 
field are independent of each other and (9) gives the Pohlhausen 
expression (reference 3) for the temperature field, which represents 
the solution for the velocity field at Pr = 1( V = a) . When the 
velocity field is known, the solution for the temperature by (9) is 
obtainable by simple quadrature . But the calculation of the velocity 
field runs into difficulties, at first, because in (8) the still 
unknown velocity appears on the right-hand side in the expression 
for f . The methods of solution by Piercy and Preston proceed from a 
random approximation for ai with which f and <l( I ) in (8) are 
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computed . The improved value cd obtained forms the starting point for 
the next step, etc. Figure 1 represents the several steps of this 
approximation method . The intentionally rough approximation to = 1 
over the entire boundary layer was chosen as orig inal solution; the 
corresponding first approximation^ to is given by the error integral . 
After the third approximation, the shearing stress Bhows a mere 
difference of 4.5 percent from the exact value. Instead of continuing 
the process mechanically, the final solution to be expected was 
estimated from the variation of the previously computed approximations 
and utilized as basis for. the subsequent step; the solution ca 

contained but a -^-percent error in shearing stress . 

With this method of solution, the improvement effected by each 
step can be est ima ted according to order of magnitude. The 
equations (8) and (9) are identical for constant density and Pr = 1. 
Assuming that the approximate solution for oj was such that for each 
individual value cu the corresponding | coordinate differed by a 
constant factor jc from the | coordinate of the exact solution, the 
effect of factor £ is then obviously just as great as that of 
quantity Pr for the temperature field. Pohlhausen found, on the 
basis of his numerical calculations, that the heat- transfer coefficient 

is proportional to thus the shearing BtresB at the wall is 

afflicted at each new step by an error of only about one-third of the 
error of the preceding step. 

For variable density, the discussed solution steps of "mathe-^ 
matical nature" can be combined with the steps of "physical nature : 

Step 1: as starting point the known solutions for constant 

density are assumed: 

(a) The Blasius solution (reference 2.) for the velocity 

profile 

(b) Pohlhausen’ s method for the .temperature field 

Step 2: 

(a) Calculation of velocity profile by (8), the temperature 

variation being based on the density of the tempera- 
ture profile according to step 1(b) 

(b) Calculation of temperature field by (9) with the velocity 

profile according to step 2(a); relation of density to 
temperature as in step 2(a) 

lit took a subsidiary worker 10 hours to reach the final solution of 
the velocity field in figure 1. 
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The process is repeated till the final solution is sufficiently 
exact, usually requiring three to four steps. 

A fey general remarks about the Influence of the temperature 
variability of the physical properties .- The flow with constant 
physical properties can be regarded as first approximation, and the 
problem is then to ascertain the differences which are produced by 
variable physical properties. The. quality of this approximation 
depends, of course^ on the temperature assumed for the physical 
properties at the isothermal flow. Choosing the wall temperature or 
the temperature at the edge of the boundary layer as reference tempera- 
ture for the isothermal flow so results on the basis of physical 
point of view £,s well, as on the basis of the equations that an increase 
of the viscosity or density inside the boundary layer is accompanied by 
an increase in the resistance; similarly, an increase in heat con- 
ductivity and density effects a greater heat transfer. But the 
magnitude of the effect of variability of the separate physical prop- 
erties is contingent upon the ratio of the boundary-layer thickness of 
the temperature and velocity field. (The ratio of both is proportional 
according to Pohlhausen.) 


This is illustrated by the following case, which is, at the same 
time, of practical importance . The temperature boundary layer is 
assumed very small compared to the flow boundary layer; consequently, 
the variation of the physical properties within the thermal boundary 
layer can be disregarded for the shearing stress and the latter 
computed as if the temperature at the edge of the boundary layer 
reaches to the wall. The same holds for the velocity profile, with the 
exception of a small area within the thermal boundary layer, where the 
velocity profile by the viscosity variation is deformed correspondingly. 
But for the temperature profile' this area is exactly decisive . 


From the equality for the shearing stresses the velocity gradients 
at the walls are: 



for Pr 


( 10 ) 


the subscript il denotes the "isothermal" flow with the physical prop- 
erties at temperature Tj_. The variability of density is noneffective for 

the field of flow, in this instance. It can be mathematically derived from 
the formulas (8) and (9) • The ratios for the temperature field are 
discussed in the next chapter by means of the two examples. 
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IH . ITiOW AND TEMPERATURE FIELD FOR VISCOUS FLUIDS 


In accordance vith the physical properties of viscous fluids, the 
velocity and temperature field were computed on the assumption that 
only the viscosity should change with the temperature by the following 
formula 


J L 


T k + T c f 

lT + T_ 


(ID 


where h and T c are constants, chosen so as to reproduce the tem- 
perature variation as closely as possible . The subscript k is to 
he 0 or 1, depending upon the choice of the physical properties 
in the dimensionless £ . The choice was h = 3 (viscous lubricating 
oil) and the two cases of a heated and cooled plate computed with 



and 8 and Pr Q 


12.5 and 100; it thus concerned identically 


great temperature differences of the same fluid, Bince Pr Q is 

for the present formed with the physical properties at wall temperature . 


Choosing T, 


as reference temperature gives by (ll) 
(p = -LL- 


in - o 


( 12 ) 


+ 1 


The result of the calculation by the iteration method of the preceding 
section is seen in figures 2 and 3« In both graphs, the dimensionless 
wall distances | Q and |j_, formed vith and \i , are plotted to 


the scale 1:\^3 and \^5:1, respectively, so that the actual wall 
distance y is the same for both abscissas . Besides the solution co, 
which took three steps to compute, the isothermal velocity profiles 
(oi) and (o>)_ at constant density at temperature T and T n 

O 1 O -L 


shovn plotted against the dimensionless coordinates 
In the subsequent compilation 


are 

2 


lo and 

^ „ and a denote the shearing stress 

at the wall and the heat-transfer coefficient = > 


For the "isothermal" temperature profiles (0) 

o 

Prandtl numbers at temperatures T Q and T^ must be 
example, in figure 2: Pt q = 12.5 and Pr_^ = 100. 


and (0)_^, the 
inserted. For 
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( T \ and (a) are the corresponding values in Isothermal flow with 
°)l 0 

the viscosity at temperature T^ and T Q . 


Table I 



T o 

a 

Pr o 

/ dm \ 

( 

'd B \ 

(T ° } ; 

(®) 0 

Jo 

Heated wall 

0.841 

1.20 

12.5 

1.58 

1.84 

Cooled wall 

I .08 

0 .98 

100 

0.255 

3.01 


Although in both cases the thermal boundary- layer thickness is far 
from small compared to the flow boundary layer, the shearing stress 
can still be computed satisfactorily by the isothermal formula with the, 
viscosity of the wall temperature. The assumptions to equation (10) 
are thus shown for Pr > 10 . 


The conditions are more complicated for the heat-transfer 
coefficient; from ( 9 ), it follows that - the heat-transfer coefficient 

a is proportional to cr(Pr) 

a is, with high accuracy, assumed as 0.664 ^/pr. Bearing in mind that 
Pr = it follows that the heat-transfer coefficient is inversely 

proportional to the sixth root of the viscosity. Since all physical 
properties except the viscosity have been assumed constant, there 
results, when it is referred once to the wall temperature 7 the other 
time to the temperature at the edge of the boundary layer 

(13) 

Vo) 

A comparison with the foregoing -tabulation indicates that (a)^ 

supplies a poorer approximation for the heat-transfer coefficient 
than (a) Q J this is readily explained by the variation of the velocity 
profile (figures 2 and 3).- It is to be expected that the conditions 
are similar at higher Prandtl numbers . 


QL 

vx 


wherein, according to Pohlhausen, 
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Another reference point for the heat- transfer coefficient is found 
in the velocity gradient at the wall; hy (10) and allowing for (3) ? 
there follows 



0 *) 


the subscripts io and il denoting isothermal flow at tempera- 
ture T 0 and T 1 « These relations are confirmed in figures 2 and 3. 


With these formulas, limits can he indicated for the heat- transfer 
coefficients (figures 2 and 3). One is given according to (13) 
by ( a )l> because the velocity profile (m)-^ yields at all points 

higher velocities at cooled and lower velocities at heated wall. The 
other limit is given hy a velocity profile of isothermal form, where 
the abscissa scale is so modified that its gradient at the wall agrees 
with the actual velocity distribution. From the remark about the 
convergence of the method of solution in II, it follows then that the 
heat-transfer coefficient is proportional to the third root of the 
velocity gradient at the wall; for this extreme value, the second 


equation of (14) gives: WEk (a) Q . Summed up, the limits of the 

^*0 

heat- transfer coefficients, *by a change in viscosity, are 



$ 





(15) 


the upper signs applying to heated, the lower to cooled wall. 

Hence, the following approximate rule for viscous fluids (Pr>10): 
For computing the resistance, the physical properties are referred to 
the temperature at the' edge of the boundary layer; for heat transfer, 
to the wall temperature . 


IT. FLOW AND TIMPERATURE FIELD AT Pr = 0 .7 (AIR) 

FOR TEMPERATURE VARIABILITY OF EVERY PHYSICAL IROPERTY 

In the -50° to ltO° temperature range the physical properties 
of the air can be represented by the following formulas 



10 


NACA 'EM 1275 


nO.780 


p = K^T 


U = KjT 

T = temperature In absolute degrees . With 


X = E T 
3 


, 0.821 


4 ■ T ° - T 1 


there results 


• q> = jl + fl(l - 0 )| 


O .78 


and similar expressions for Tjr and X . 


and ~ 


showed 


The calculations for a heated plate and tS = -^ 

only moderate differences in the velocity and temperature field from 
the form for isothermal flow (Table 2) . For the investigation of the 
conditions at higher temperature differences, the case T^ = 20 

and T 0 = 620° C was computed. The velocity and temperature fields 
already exhibit, according to figure 4, appreciable differences from 
the form for constant physical properties; -and are formed 

with the physical properties at temperatures T Q and T^, respectively. 

This results in a substantial thickening of the boundary layer for both 
fields; nevertheless, wall shearing stress and heat-transfer coefficient 
indicate only minor departures from the values for constant physical 
properties . 


Table 2 


Heating 

®L 


T o 

( T o)o 

T o 

ror 

a 

?<3o 

a 

(a-)l 

4,1 

0.575 

0 .h 90 

1.02 

1.00 

1.01 

1.00 


•514 

.420 

1.05 

1.00 

1.02 

1.00 

T q - 620° C 
T-jj = 20° C 

.286 

.235 

1.11 

•93 

1.03 

.96 
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Tile explanation for it is that in air the growth of the viscosity with 
the temperature acts in the sense of a resistance increase, the drop 
in density in the sense of a resistance decrease, and "both effects 
practically cancel one another at Pr = 0.7, where thermal and flow 
boundary-layer thickness are about equally great. The conditions for 
the temperature field are almost identical, because the heat con- 
ductivity and the viscosity are similarly affected by the temperature. 


The frictional heat can be allowed for in similar mann er; 

^2 

equation (2) contains then an additive term p 


side, and the solution reads 



on the right-hand 


0 


A(g) 

BCi) 

H(5) 


P h 

J 0 X 


-E(!) 


d| 


_ o Fr k M 'e 
(Tl - T 0 ) 

1 


r I fr 6 qn/ 2 e R C£) d6 

Jo X \Jo 


= Pr, 


rmi 

Jn X 


de 


e “E( | ) d| 


>> (16) 


Zd = 
e 


2c p 


The iteration method can be applied again, although a little more 
paper work is involved. For constant physical properties, equation (l6) 
reduces to Eckert's solution (reference 8). The thermometer problem 
(vanishing temperature gradient at the wall) can also be solved by 
suitable variation of the integration, constant. In view of Crocco's 
calculations for a gas with Pr = O. 725 , it was decided not to cal- 
culate any model problems by the new method. 
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V . APPLICATION TO A BZETOSION PROBLEM 


The concentration field for the problem of diffusion at the flat 
plate can "be calculated in the same manner as the temperature field.3>^- 
The differential equation reads “ 


2 

u^£ + v— = k^— 2. 
ox by Sy2 


(IT) 


where k is the diffusion factor and c the concentration which is 
defined as quantity of gas or vapor per unit volume . The physical 
properties are regarded as constant, hut it is also taken into account 
that for greater concentrations the velocity v at the wall no longer 
disappears, as already pointed out hy Nusselt (reference 7)* 'When 
fluid from a wall is vaporized, say hy a gas such as gas flowing along 
a wetted wall, substance passes continuously into the flow. 

Hence v(0)>0 at the wall. When, on the other hand, vapor condenses 
at the wall or when air containing ammonia, for example, passes over 
blotting paper impregnated with hydrochloric acid, it results 
in v(0)<0 v The houndary' conditions for v are according to the 
equations (lOO) and (101) of reference (7)s 


- 1 = T (0) ( l8 ) 
co\oy J Q x 

Po 

where c is the concentration of the gas or vapor, for which the wall 
is permeable, c 0 the concentration at the wall, p Q the corresponding 

partial pressure, and p the total pressure . _ 

Introduction of the flow velocity U and the dimensionless | 
results in —•••-- 

/ v\ = _ 1 Js. C 1 " c o /dc ^ .AT c = 0 " °o (19) 

\U/ o 2 Co _ 1 j u \dg/ 0 Vvx __ c 1 - c 0 

^Eckert reported a solution of this problem at the 19^3 meeting 
of the TDI committee for heat research in Bayreuth, where - an approxima- 
tion method similar to Pohlhausen's method for the flow boundary layer 
was used . 

It ii . 

Damkohler s estimate for the present problem was published in 
Z. fur Slektrochemie, 19^2, p. 178. _ ... 
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where c D and c^ are the concentration at the wall and at the edge 
of the boundary layer. Similarly to (5), it gives 



and similarly to (9)» the solution for the concentration field 




' if* 


M)d| 


<4 


f(l) = 2 £ d| 

d n U 


> ( 21 ) 


where !. is a quantity analogous to the Prandtl number. To obtain the 
k 

velocity cn, simply put 3d » 1 in equation (21). 

k 


The concentration gradient at the wall is contained in M; but (21) 
can, in the first instance, be solved for any M values and the quantity 

computed with the aid of the value obtained from 


calculated M and IT values are represented in figures 5 and 6, the 
concentration gradient at the wall, in figure 7 • M>0 denotes evapora- 

tion at the plate; M<0, condensation and absorption at the plate; the 

value 0 .6 chosen for the quantity 3d is applicable in good approxima- 

k 5 

tlon for the diffusion of water and ammonia in air. Strictly speaking, 
for the specified higher concentrations, the density and viscosity of 
the two fluids are dependent on the concentration; and the diffusion 
factor, on the temperature. Cases of that kind can be calculated 
with the aid of the described method. If the diffusion is 


1 • The velocity and concentration fields for the 
o 



pp. 


''According to Ten Bosch: 
189 and 257. 


it ti 

Die Warmeubertragung, Berlin 193° , 
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accompanied by a heat transfer, the solution for the concentration field 
can equally he applied to the temperature field with good approximation* 
In the same vay, the heat transfer can he derived from the solution for 
the concentration field, when air is exhausted or blown at the plate 
with transverse velocities at the wall corresponding to. equation (20). 


CONCLUDING NOTE 


After completion of the calculations the writer received knowledge 
of a report hy Schlichting and Bus smarm (reference 11) about the 
velocity profile at the flat plate for exhaustion where the transverse 
velocity at the wall was expressed hy 



Between the present value M and C the following relation exists 
(see also (19) and (20)). 

C = -M 

The present velocity distributions agree to about 1 percent with those 
calculated by Schlichting (by a different method), with exception of 
M = 1, where the writer plainly chose too few approximation steps and 
the differences are therefore a little greater. For the present 
calculation three steps were usually sufficient. 


Translated by J. Yanier 
National Advisory Committee 
for Aeronautics 
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1.0 2.0 3.0 


Figure 1.- The several approximations for computing the velocity- 
distribution at the flat plate by the method of Preston and Piercy 
(constant physical quantities). 



Figure 2.- Velocity and temperature distribution at a heated plate 
for variable viscosity. Viscosity exponent b = 3. (w) 0 , (e) 0 , 

and (u)^, (e)^ isothermal velocity and temperature distributions, 

v Q and kinematic viscosity at wall temperature T Q and 

temperature T]_ at the edge of the boundary layer. 
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Figure 3.- Velocity and temperature distribution at the cooled plate. 
■ Viscosity exponent b = 3. (u) 0 , (e) Q , and (w)^, (e)^ isothermal 

velocity and temperature distributions, v Q and v-^ kinematic 
viscosity at wall temperature T 0 and temperature at 
the edge of the boundary layer. 



Figure 4.- Velocity and temperature distribution at a heated plate 
for Pr = 0.7 (air); all physical quantities constant with 
temperature. 
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Figure 5.- Velocity field at diffusion with higher concentrations, where 
finite transverse velocities occur at the wall (see text for equations (18) 
to (20)). 



Figure 6.- Concentration distribution to figure 5. 
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Figure 7.- Concentration gradient at the wall and quantity M plotted 
against N (see text for equations (19) and (20)). 
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